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Abstract 

The Df — > K ± K t tt ± absolute branching fraction is measured using e + e~ — > D* ± D^(2536) 
events collected by the Belle detector at the KEKB e + e~ asymmetric energy collider. Using the 
ratio of yields when either the D s \ or D* is fully reconstructed, we find B(Djr — ► X^-K^tt^) = 
(4.0 ± OA(stat) ± 0A(sys))%. 

PACS numbers: 14.40.Lb, 13.66.Bc, 13.25.Ft 
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Knowledge of meson 1 absolute branching fractions is important for normalization of 
many decays involving a D t in a final state. The poor accuracy of the branching fraction 
B(D+ -> K+K--K+) = (5.2 ± 0.9)% [l| has been a systematic limitation for some pre- 
cise measurements. In particular, the recent study of the CP violation in B° — > D^ ± n T 
decays is restricted by the knowledge of the ratio of two amplitudes that determine the 
CP-asymmetry 2], ^]. The amplitude B° — > D^ + ti^ can be calculated from the branching 
fraction of B° — > D^ + n~ decays assuming factorization. On the other hand, the factoriza- 
tion hypothesis can be tested by measuring the ratio of B° -> D^-ji + and B° -> D^~D+ 
decays. Both B(B° -> D^+n') and B(B° -> D^~D+) measurements can be improved 
with better accuracy in absolute branching fractions. 

Recently, the absolute branching fraction of — > 0(— > K + K~)n + was measured by the 
BaBar collaboration, which used partial and full reconstruction of B — > £)(*)£)(*)+ decays 4]. 
Another result obtained from a y^-scan above D+D~ threshold was presented by the CLEO- 
c collaboration [||. 

In this paper we report on a measurement of the — > K + K~7T + branching fraction 
using two body e + e~-continuum annihilation into a -D* + _D,Ti(2536) final state. The analysis 
is based on 552.3 ftT 1 of data at the T(4S) resonance and nearby continuum, collected with 

n n 

the Belle detector [6( at the KEKB asymmetric energy storage ring [7j. 
I. METHOD 

We use the partial reconstruction of the process e + e~ — > D* + D~ x . In this analysis 4- 
momentum conservation allows us to infer the 4-momentum of the undetected part. The 
method used was described in Ref. 8| and applied to the measurement of the e + e" — > 
£)(*)+_£)(*)- cross sections. 

Here we reconstruct the process e + e~ — > D^D^ using two different tagging procedures. 
The first one (denoted as the D~ x tag) includes the full reconstruction of the D~ x meson via 
— > _D*.K" decay and observation of the photon from Z}* + — > .D+7, while the Z}+ is not 
reconstructed. The measured signal yield with the Dj± tag is proportional to the branching 
fractions of the reconstructed D* modes. In the second procedure (denoted as the D* + tag) 



1 Charge conjugation is implied through the paper. 
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we require full reconstruction of D* + through D* s + — > -D+ 7 and observation of the kaon 
from — > but the D* is not reconstructed. Since the Df meson is reconstructed 

in the channel of interest, Df — > K + K~ir + , the signal yield measured with the D* + tag 
is proportional to this branching fraction. The (efficiency-corrected) ratio of the two 
measured signal yields is equal to the ratio of well-known D* branching fractions and the 
branching fraction of the Df. 

In order to calculate reconstruction efficiencies and optimize event selection criteria, 
Monte Carlo (MC) signal events are generated and simulated using a GEANT-based full 
simulator, including initial state radiation (ISR), and assuming no form-factors for D* + and 
D~ x mesons. 

To identify the signal we study the mass recoiling against the reconstructed particle (or 
combination of particles) denoted as X. This recoil mass is defined as: 

M rccoil (X) = ^f(E CM -E x ) 2 -P x , (1) 

where E x and Px are the center-of-mass (CM) energy and momentum of X, respectively; 
E C m is the CM beam energy. We expect a peak in the M recoi i distribution at the nominal 
mass of the recoil particle. 

The resolution in M recoi i is ~ 50MeV/c 2 according to the MC, which is not sufficient to 
separate different final states, e.g. D^D^, D^D^ and non-resonant D^DK. To disentangle 
the contribution of these final states we use another kinematic variable, the recoil mass 
difference AM rccoi i: 

AM recoa ( J D a - l7 ) ee M iecoil (D~) - M^D;^) , (2) 
AM rccoil (D; + ^) ee M recoil (L>:+) - M rccoil (D* s + K). (3) 

In the tag procedure the signal events make a narrow peak in the AM reco ii(-D7i7) dis- 
tribution at the nominal D* + — D+ mass difference with a resolution of a ~ 5MeV/c 2 , 
dominated by the 7 energy resolution according to the MC. The AM recoi i(D* + i^) spectrum 
peaks at the — D* mass difference with a resolution of a ~ 2 MeV/c 2 . 

II. SELECTION 

All primary charged tracks are required to be consistent with originating from the in- 
teraction point. Charged kaon candidates are identified using information from dE/dx 
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measurements in the central drift chamber, aerogel Cherenkov counters and time-of-flight 
system. No identification requirements are applied for pion candidates. Kg candidates are 
reconstructed from 7r + 7r~ pairs with an invariant mass within 15MeV/c 2 of the nominal K s 
mass. The distance between the two pion tracks at the K s vertex is required to be smaller 
than 1 cm, the flight distance in the plane perpendicular to the beam from the interaction 
point is required to be greater than 0.1 cm and the angle between the K s momentum direc- 
tion and decay path in this plane should be smaller than 0.01 rad. Photons are reconstructed 
in the electromagnetic calorimeter as showers with energies above 50 MeV that are not asso- 
ciated with charged tracks. n° candidates are reconstructed by combining pairs of photons 
with invariant masses within 15MeV/c 2 of the nominal 7T° mass. 

D° candidates are reconstructed using five decay modes: K~tt + , K~K + , K~7t~7t + tt + , 
K s n + n" and K^n + n°. D + candidates are reconstructed using the K~it + it + channel. Df 
candidates are reconstructed using the K + K~n + decay mode. A ±15MeV/c 2 mass window 
is used for all D modes (approximately 2.5a in each case) except for the D° — > K~it~it + it + 
where ±10MeV/c 2 is applied and D° -> A^Tr+Tr with ±20MeV/c 2 mass window. All £>+, 
D + and D° candidates are subjected to a mass and vertex constrained fit to improve their 
momenta and thus the recoil mass resolution. 

D* candidates are selected via D* + — > D°ir + and D*° — > D°tt° decay modes with 
±2MeV/c 2 mass window. D* + candidates are reconstructed via channel with 

±10MeV/c 2 mass window. is reconstructed in D*°K~ and D*~K S decay modes. 

In the case of multiple candidates, the candidate with the minimum value of xfot (Xtot = 
Xm(d) + Xm(d*) for D 7i ta S and Xtot = Xm (d +) for D * s + tag, respectively) is chosen. Each x 2 
is defined as the square of the ratio of the deviation of the measured mass from the nominal 
value and the corresponding resolution. 

III. RECONSTRUCTION 
A. tag 

As the ratio of D~ x — > D*°K~ and — > D*~K S branching fractions is unknown, we 
perform the analysis for these two channels separately and average the results. 

The D*°K~ and D*~K S mass spectra with a preselection requirement on M rccoi \(D* K) < 
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3 GeV/c 2 are shown in Fig. [TJ Clean peaks at the nominal D~± mass without substantial 
combinatorial background are evident in both distributions. We fit these distributions with 
a sum of two Gaussians with a common mean and the background parametrized with a 
threshold function: 

Bg{m) = (A ■ m + B) ^Jm - Mgf G - M| DG . (4) 

The fit yields a Gaussian central value of M = 2535.1 ± 0.4MeV/c 2 for D*°R- and M = 
2535. 0±0. 3 MeV/c 2 for D*~K$ channels, respectively, which is in good agreement with both 
the PDG and MC values. 




FIG. 1: The distribution of invariant mass a) D*°K and b) D* K$. The solid curves represent 
the fit described in the text. 

The signal and sideband regions for D~ x are defined as \M{D~i) — 2.5354| < 0.005 GeV/c 2 
and 0.010 GeV/c 2 < |M(Z?7i) — 2.5354| < 0.015 GeV/c 2 , respectively. The M recoil (.D;i) spec- 
tra are shown in Fig. [2j The significant excess in the region from 1.9 to 2.2 GeV/c 2 includes 
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e + e~ — > DfD~ x and e + e" — > D* s + D~ x processes. The small excess around 2.3GeV/c 2 could 
be interpreted as a contribution from the process e + e~ — > D^jD^. The excess around 
2.5GeV/c 2 corresponds to the excited D** + = D^, Df 2 states. 
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FIG. 2: a) M recoa (D7 1 -» D*°R-) distribution; b) M recoi i(D^ -> Z>*~.fifg) distribution. Points 
with error bars correspond to the events in the signal region. The sideband contribution is 
shown as hatched histogram. The signal windows are indicated by the vertical lines. 

We define the signal M reco ii(-DJi) window as \M recoil (D^) - 2.1 121 1 < 0.100 GeV/c 2 
(Fig. [2]). As the recoil mass resolution of ~ 50MeV/c 2 is not sufficient to resolve 
e + e~ — ► D+D~ Y and e + e~ — > D* s + D~ x processes, to disentangle them we reconstruct the 
7 from D* + — > D+j decays. The reconstructed D~ ± candidates from the signal M reco ii(-DJi) 
region are combined with each photon in the event to calculate the recoil mass difference 
AM recoi i(Lr l7 ). 

To obtain the signal yield and disentangle the signal contribution from the non-resonant 
D* + D*K process, we perform fits to D*K mass distributions in bins of AM recoil (D7i7)- The 
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AM reco u(-DJi7) bin width is chosen to be 4MeV/c 2 . The signal function is the sum of two 
Gaussians with the mass and width fixed to the corresponding MC values. The background 
is described with a threshold function (Eq. H]). The signal yield for each bin of the recoil 
mass difference for the D~ x — > D*°K~ and — > D*~K$ channels are plotted in Fig. [31 
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FIG. 3: The signal yields in bins of AM reco u(D sl 7): a) for the D sl -> D*°K~ channel; b) for the 
D~i —* D*~Kg channel. The solid curves represent the fits described in the text. 

The method used to extract the signal yield guarantees that the signal peak includes 
only those events that contain both a D~ x meson and a D* + decaying to D s <y. Thus, the 
only possible peaking background sources are final states with extra pion(s) (e.g. e + e~ — > 
D^D^n ). The process with one extra tt° is suppressed due to isospin conservation. The 
possible contribution of e + e~ — > -D7i-Djj(2460) followed by D (2460) — > D* + n° decay, where 
the isospin violation occurs in the D s j decay, is suppressed by the recoil mass cut at the 
4cr level. The processes with two pions are also suppressed at the ~ 3.5<r level. Therefore, 
we conclude that the remaining background is negligibly small; the extracted yields are 
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dominated by the signal. 

The distributions in Fig. [3] are fitted to the sum of a signal Gaussian (with the mass and 
width fixed to the MC values) and a linear background function. The fits yield 184 ± 19 
signal events for — > D*°K~ and 105 ± 14 events for 

reconstruction efficiency the same procedure is applied to the MC signal sample. 



D* Kg. To determine the 



B. D*+ tag 

In the D* + tag, the D* + of the process e + e~ — > D* + D^ is fully reconstructed. The 
.D+7 mass distribution is shown in Fig. H] and fitted with the sum of a signal Gaussian and 
a linear background function. The obtained central value of the Gaussian M = 2113.5 ± 
0.4MeV/c 2 is somewhat larger than the PDG value, but in good agreement with the MC 
expectation, which shows a similar shift. The signal region is chosen to be \M(D^j) — 
2.1121| < 0.010 GeV/c 2 while the sidebands are defined as 0.015 GeV/c 2 < \M(D+i) - 
2.1121| < 0.025 GeV/c 2 . 
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FIG. 4: The Df"f invariant mass distribution. The solid curve represents the fit described in the 
text. 



The M recoi i(D* + ) spectrum for the signal region is shown in Fig. [51 A wide enhance- 
ment around 2.1 GeV/c 2 is formed by overlapping D~ and D*~ signals. The signal process 
e + e~ — > £)* + £)7i contributes to the bump around 2.55GeV/c 2 together with the process 
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e + e — > D* + D s2 . We define the signal region by the requirement \M recoi \(D* + ) — 2.5354| < 
0.100 GeV/c 2 (~ 2a). 



2000 




FIG. 5: The M reco n(D* + ) distribution. The events in the D* + signal region are shown as points 
with error bars. The D* + sidebands are superimposed as hatched histogram. The M reco ii(L>* + ) 
signal window is indicated by vertical lines. 

To extract the e + e~ — > D* + D~i signal yield for D* + tag, we perform a fit to the -D^7 
mass distributions in bins of AM Tecoi \(D* + K) . The signal function is a Gaussian with its 
parameters fixed to the corresponding MC values. The background is parametrized by a 
linear function. The yields returned by the fits for D~y — > D*°K~ and Dj x — > D*~Kg 
channels are plotted in Fig. [6j We fit the resulting distribution with the sum of a signal 
Gaussian (with fixed mass and width) and a threshold function (Eq. HJ). The fits yield 
267 ± 18 signal events for the D*°K~ channel and 45 ± 8 events for the D*~Kg channel. 

IV. SYSTEMATIC ERRORS 

The systematic uncertainties are summarized in Table [B 

Different production and helicity angular distributions result in different ratios of re- 
construction efficiencies for the two tagging procedures. A potential source of significant 
systematic error are the uncertainties, associated with reconstruction efficiencies of the kaon 
from the D~ x decay and the photon from the D* + , which strongly depend on the mother 
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FIG. 6: The e+e~ D^+D^ yields for D*+ tag in bins of AM Tecoil (D* + K): a) for the D~ x -> 
D*°K~ channel; b) for the -» D*~ K% channel. The solid curves show the fits described in the 
text. 

particle polarizations. As the kaon and photon are, however, reconstructed in both tags, the 
systematic uncertainties related to the reconstruction of these particles cancel out. The re- 
maining systematic uncertainty in the efficiency ratio from the production angle # pro d (polar 
angle of the momentum of D~ x meson) is estimated from MC simulation by comparing the 
regions | cos(6 , pro d)| > 0.5 and | cos(0 pro d)| < 0.5 and found to be small (0.4%). 

The systematic error arising from unknown D* polarizations contributes only to the 
uncertainty in the efficiency of D~ t tag and does not cancel out in the efficiency ratio. This 
contribution is estimated from MC assuming different D* polarizations. 

The efficiency of the requirement on recoil mass depends on the fraction of events with 
an ISR photon as well as on the recoil mass window cut. The error from the ISR photon 
is conservatively estimated by changing the fraction of events containing ISR photons in 
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TABLE I: Contributions to the systematic error. 



Source 


Estimated error, [%] 


Production angle 


0.4 


D* polarization 


3.8 


ISR 


2.3 


Recoil mass cut 


4.8 


Fit 


4.3 


Tracking/7r° reconstr. efficiency 


3.0 




5.2 


MC statistics 


2.3 


Total 


10.1 



the whole MC sample by 30%. The resulting relative shift in B(Df) is found to be 2.3%. 
The error from the recoil mass cut is estimated by varying the recoil mass window width 
in the range of 80 - 120MeV/c 2 . The difference in the calculated value of B(D+) (4.8%) is 
conservatively considered as the systematic error. 

The fit systematics is estimated assuming different shapes for the signal (single and 
double Gaussians) and for the background (different order polynomials). The difference 
in the resulting branching fraction is found to be 4.3% and is treated as the systematic 
uncertainty in the fit. 

Other contributions come from the uncertainties in the track/V reconstruction efficien- 
cies, branching fractions and MC statistics. 

All of the systematic contributions described above are added in quadrature to give 10.1%. 



V. RESULTS 



Using the measured signal yields N(D* + ) and N(D~ 1 ) with D* + and D~ Y tags, respec- 
tively, and taking into account the efficiency ratio ( J p * + ) = 0-50 for D*°K~ and ^§Srj — 1-65 
for D*~Kg channels, respectively), we obtain the Df absolute branching fraction: 

B w- m+ >=wlly' B(I7W) ' (5) 
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where B(D^) is the product of D* branching fraction and those of sub-decays. Separate 
calculations for the D* R- and D*~K° S channels yield B(D+ -> K + K-n + ) of (4.01 ± 
0A7(stat))% and (3.84 ± 0.83(stat))%, respectively. 

VI. SUMMARY 

In summary, we report a measurement of the —>■ K + K~n + branching fraction using 
a new method of double tag partial reconstruction. The branching fraction is measured 
separately in two channels e+e~ -> D* s + D^(-> D*°K~) and e+e" -> D* + D^(^ D*~K° S ). 
The average value is B(D+ -> K+R-7T+) = (4.0 ± 0.4(stat) ± 0.4(sys))%. 
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